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The generation of high harmonics in laser-irradiated light homonuclear diatomics (-^2^, N2 
and O2) compared to that in atomic counterparts (of nearly identical binding energy) is stud- 

J2li| ied within the velocity gauge version of conventional strong-field approximation. The applied 

strong-field approach (alternatively developed earlier to incorporate rescattering effects be- 
yond the conventional saddle-point approximation) is currently extended to molecular case by 

j^ . means of supplement the standard linear combination of atomic orbitals and molecular orbitals 

c/j ! method. The associated model proved to adequately reproduce a general shape and detailed 

structure of molecular harmonic spectra, which demonstrate a number of remarkable distinctive 

►l^| differences from respective atomic spectra calculated under the same laser pulses. The revealed 

J^ ■ differences are found to be strongly dependent on internuclear separation and also very sensitive 

"■ to the orbital and bonding symmetry of contributing molecular valence shell. In particular, the 

model correctly predicts the behavior of high-frequency plateau (both for its extent and even 
details of structure) in molecular harmonic spectra at small (nearly equilibrium) and large in- 
ternuclear separations. In addition, for some group of harmonics, the harmonic emission rates 
were ascertained to dominate by contribution from inner molecular shells of higher binding 

Tj- . energy and different orbital symmetry compared to the outermost molecular orbital normally 
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predominantly contributing. 
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1. INTRODUCTION AND BACKGROUND MOTIVATION. 

Apart from pure fundamental interest, the strong-field phenomenon of high-order 
harmonic generation (HHG) in laser-irradiated species is currently believed to be a highly 
promising for creation compact (table-top) sources of ultra-short (sub-femtosecond) pulses 
of powerful coherent XUV radiation [1, 2]. Being accessible, such a high-frequency radia- 
tion might be used for X-ray holography or monitoring of biological processes in vivo [2] 
as well as tomographic imaging of molecular orbitals [3] with an unprecedented brightness 
and spatial and temporal resolution. For the reasons above, the HHG process became 
a subject of equally great experimental and theoretical interest in the last decade (see 
also e.g., [4-6], for review and references cited therein). Compared to atoms, the laser- 
irradiated molecules as sources of high harmonics [7-14] offer promising prospects as to 
raising the up to now disappointingly low conversion efficiency for harmonic generation 
[7-12] and/or harmonic cutoff frequency [10, 13] due to an advanced possibility of con- 
trolling and manipulating the additional external degrees of freedom in molecules. For 
example, the internuclear distance Rq of simplest diatomic molecules, which determines 
differently delocalized initial electronic states, can be used to enhance the harmonic con- 
version efficiency exceeding the atomic one due to transient enhancement of harmonic 
emission in expanding molecules [12]. Particularly, the calculated harmonic conversion 
efficiency was found to be noticeably enhanced (by several orders of magnitude) with 
increasing of Rq over equilibrium internuclear separation [10, 12, 14]. Diatomic molecules 
may also produce high harmonics of much higher cutoff frequency due to a much longer 
high-frequency plateau extended up to 8Up of maximum length [10] versus the respective 
maximum value 3.17Up ascertained for atomic harmonic spectra (here Upis the pondero- 
motive energy of field-induced oscillating motion of an electron in the continuum; the 
atomic system of units is used unless stated otherwise). 

Because of the electron wavelength is comparable to the size of the molecular wave- 
function, the intensity of emission at a particular X-ray photon energy is also expected 
to be sensitive to the wavelength of the electron wave and the shape of the molecular 
wavefunction and/or orientation of internuclear molecular axis with respect to incident 
laser field polarization [8-11]. In addition, simultaneous observation of both ion yields and 
harmonic signals [15] under the same conditions was shown to serve as an effective tool 
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to probe molecular structure in an instant [16] and also make it possible to use quantum 
phenomena in HHG associated with molecular symmetries to disentangle the contribu- 
tions from the ionization and recombination processes [17]. In this context, homonuclear 
(or homopolar) diatomics seem to be also of special interest due to identical atomic 
nuclei give rise to a possibility of various laser field-induced and orientation-dependent 
intramolecular interference effects [15] that may additionally affect the molecular HHG 
process. The clear signs of intramolecular interference-related effects due to two-centered 
nature of diatomic molecule were numerically revealed [11] in high-harmonic spectra of 
model [ID or 2D) H2 and H2 systems as a broad and pronounced local minimum of 
the position strongly dependent on the internuclear separation and spatial orientation of 
diatomic molecule. The origin of the minimum was attributed to the destructive interfer- 
ence of photoelectron emission from separate atomic centers [11], so that a transparent 
interpretation can be given in terms of interference between two radiating point sources 
located at the positions of the nuclei. The latter pivotal idea (obviously proposed first 
in Ref. [18]) about intramolecular interference of atomic photoionization amplitudes was 
successfully exploited in [19] for description of molecular ionization as compared to ion- 
ization of respective atomic counterparts having nearly identical binding energies. Based 
on the velocity gauge formulation of conventional strong-field approximation (SFA) (see, 
e.g. [20] and associated Keldysh-Faisal-Reiss theories [21]) with the atomic ionization 
rates modified by the interference from the atomic centers, the molecular ionization was 
predicted to be enhanced (or, at least, not suppressed) in N2 (versus Ar) and/or sup- 
pressed in O2 (versus Xe ). The mentioned suppression in O2 was interpreted in terms of 
the strong-field MO-SFA model [19] as intimately related to intramolecular interference, 
which is to be always destructive for low-energy photoelectrons emitted from the ItTq 
highest- occupied molecular orbital (HOMO) of antibonding symmetry and predominantly 
contributing to respective ionization rate. Accordingly, the background physical mecha- 
nism underlying an enhanced ionization in N2 versus Ar was identified as constructive 
intramolecular interference for low-energy photoelectrons emitted from the Sa^ HOMO 
in N2 of bonding symmetry. 

Meantime, because of the high suppression in strong-field ionization of O2 versus 
Xe, the laser-irradiated diatomic O2 is expected to produce high harmonics of notice- 
ably higher cutoif frequency as compared to Xe atom [13]. Namely, due to a suppressed 
molecular ionization in H2 and O2 (versus their atomic counterparts) and correspond- 
ingly larger value of the saturation intensity, the latter diatomics survive at a consid- 
erably higher laser intensity (at which their atomic counterparts are to be completely 
ionized) and, thus, they are still able to produce harmonics of a considerably higher cut- 
Oj^ frequency (proportional to laser intensity). The related extension of high-frequency 
plateau in harmonic spectra of O2 relative to Xe, was recently observed in experiment 
and numerically simulated [13] based on the Lewenstein model [22] supplemented by the 
molecular tunneling ionization (or the so-called MO-ADK) theory [23] to account for the 
differences in ionization of molecules compared to atoms in terms of Ammosov-Delone- 
Krainov ( ADK) model of tunneling ionization [24] . The MO-ADK model partly succeeded 
in reproducing of suppressed ionization (e.g., in O2 versus Xe compatible with experiment 
[25]), although, being applied to molecular HHG, the model seems to encounter difficulties 
in providing a transparent interpretation and satisfactory agreement with experiment in 
prediction of the cntojff frequency observed in molecular high-harmonic spectra [13] (see 
also Sec. 4 below, for details). Being also applied to ionization of other diatomic species 
(e.g., in H2 versus companion Ar atom), the ADK-based approach resulted in a dis- 



crepancy with experiment [26], at least in the strong- field domain (i.e., only where the 
tunneling theory is generally valid). Since the wavefunction of initial molecular state cor- 
responding to HOMO is always approximated by only a single one- centered atomic orbital 
within the MO-ADK model, such theory is not suitable either for adequate description 
of intramolecular interference-related phenomena arising from the two-centered nature of 
molecular valence shell in diatomics. 

Thus, there are numerous interrelated and remarkable aspects of molecular HHG 
and above-threshold ionization (ATI) as a prerequisite to the HHG, which appear to be 
also strongly interdependent and interrelated in laser-irradiated molecules too. However, 
with the possible few exceptions, the fundamental physics associated with how atoms and 
molecules respond to intense laser radiation (i.e. optical fields approaching the atomic 
unit) cannot be treated exactly due to the comparable strengths of the electron-nucleus, 
electron-electron and electron-field interaction. Although ab initio quantum calculations 
for atoms are readily available, at least within the single-active electron (SAE) approxima- 
tion, this is not the case for molecules as fully three-dimensional {3D) calculations on even 
relatively simple diatomics, such as H2, are extremely difficult [14]. Therefore, the HHG 
process in diatomic molecules is primarily treated nowadays using approximate SAE- 
based general strong-field approaches, such as conventional strong-field approximation 
(SPA) formulated in terms of "rescattering" (or recollision) picture and model zero-range 
binding potential of atomic centers [10] and/or various pure numerical procedures (see e.g. 
[27]) and simulations (e.g., [8, 11, 12]). The harmonic emission rates for molecules can be 
also calculated, in principle, based on different ab initio numerical procedures and meth- 
ods, such as direct numerical solution of time-dependent Schrodinger equation (TDSE 
method, e.g., [11]) and time-dependent density-functional theory (TD-DFT) including 
many-electron effects [28]; however, the related results seem to be very computationally 
demanding and hardly available for reliable and transparent interpretation. 

All the aspects above gave us a strong motivation to consider the harmonic genera- 
tion in laser-irradiated diatomics and verify the mentioned aspects within the currently 
proposed alternative model of molecular HHG. The model is based on fully quantum- 
mechanical SAE-based strong-field approach developed earlier [29, 30] and proved to ad- 
equately incorporate rescattering effects in atomic HHG and high-energy ATI processes. 
In the present paper the proposed strong-field approach is extended to molecular case 
of diatomic specie (i.e. substantially many-electron system) by means of supplement 
the standard linear combination of atomic orbitals (LCAO) and molecular orbitals (MO) 
method invoked for approximate reproducing the two-centered one-electron wavefunction 
of initial molecular discrete state. The LCAO-MO method proved to be surprisingly 
well working for adequate description and transparent interpretation the phenomena of 
enhanced (and/or suppressed) ionization of laser-irradiated diatomics in terms of pre- 
viously developed strong- field SFA-LCAO model [31]. This model is to be also used 
further in Sec. 3 for numerical calculation and evaluation of laser intensity for ionization 
saturation under conditions of molecular HHG experiment [13]. The related results for 
calculated harmonic emission rates (or high-harmonic conversion efficiency) and position 
of high-frequency plateau cutoff proved to be sensitive to the molecular orbital symmetry 
and internuclear separation. At last, our present consideration is not restricted to only 
one, a single (the outermost) molecular valence shell, so that other (inner) contributing 
molecular orbitals of higher binding energies and different orbital symmetry are addition- 
ally and equally well incorporated, wherever appropriate. The proposed molecular HHG 
model is shown to adequately reproduce the distinctive features of molecular harmonic 



spectra compared to respective atomic spectra as well as noticeable differences found in 
structure of molecular harmonic spectra corresponding to HOMO of bonding symmetry 
(e.g., Sag in A^2) relative to antibonding symmetry (e.g., liTg in O2). 

2. THE APPLIED STRONG-FIELD HHG MODEL: BASIC ASSUMPTIONS AND ANA- 
LYTICAL RELATIONS. 

The harmonic emission process under discussion is to be further investigated within 
framework of the strong-field fully quantum-mechanical approach developed in [29, 30] 
and applied previously to treatment of atomic HHG process. Let's outline first in brief the 
main features of related strong-field model, which is currently to be extended to molecular 
case. Within the dipole approximation (neglecting any photon momenta) the Hamiltonian 
of EM interaction in the velocity gauge (VG) form, reads as 

W{r,t) = -A{t)-p + ^A'{t) (1) 

Co ZCq 

where p= — iV is the operator of electron canonical momentum, cq ~ 137 is the light 
velocity in vacuum. This means that the vector potential A(t) of incident laser field and 
associated field strength E(t) are supposed to be independent on coordinate r, but both 
are functions of time t only. Our present consideration will be restricted entirely to the 
case of only harmonic emission processes in which the atomic or molecular specie occupies 
the same state before and after the passage of the laser pulse. Thus, the bound discrete 
state is supposed to be unperturbed by incident laser field and the respective undis- 
torted (laser-free) stationary wavefunction $0 (r, t) = $0 (r) exp {—ieot) already known in 
advance (or, at least, found to any arbitrary prescribed accuracy). 

Then, for particular case of linearly polarized laser field of the frequency uj, electric 
field strength E and unit polarization vector e: 

Al (t) = (co/cj) eE cos (cut) (2) 

the total SFA-amplitude of HHG can be expressed through the sum of partial amplitudes 
/n (^k') of emission of high-harmonic of A^-th order, frequency i7k' and polarization 
ek',A' (A = 1,2): 

/i^r^ (fikO = E fN^'''^ (^kO S (l^k' - Nuj) = (3) 

N 
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Qn (t) = tn^'Dn (t) = 2fik' E ^n^ (^) exp {-tNoot) (4) 

N 

with the time-dependent dipole moment Dq (t) induced in a laser-irradiated system by 
incident laser field, whereas V is the spatial normalization volume. The A^-th Fourier 
component Dq (cu) of the field-induced dipole moment can be derived in the following 



analytical form 

Dq i^) = E RN-m,m{^k',qm,il) = -'^'^^k^ E {Up - muj) X (5) 
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Here 

Fo (q) = {2Try'/^ J dv exp (-^q ■ r) <l>o(r) (6) 

is the Fourier transform of stationary wavefunction $o(i") of laser-exposed specie corre- 
sponding to initial bound discrete state unperturbed by incident laser field, while Bg (x; y) 
is the s-th order generalized Bessel function [20] of two real arguments as dimensionless 
parameters: 

C(p) = (E■p)/cu^ r, = U,/u = E^/(Auj^) (7) 

with ponderomotive energy Up of field-induced oscillating motion of electron driven in 

continuum. The variable q^ = J '2 {rnu — Ip — Up) in Eq.(5) denotes discrete values of 
photoelectron canonical momentum in intermediate continuum states of energies 
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M =muj-Ip-Up (8) 



corresponding to m number of absorbed laser photons beginning from the minimum one 
A^o = ^0 iv) = [{^p + Up) /uj] + 1 required for ionization (here [x] denotes an integer part 
of variable x) . 

Due to the "pole approximation" applied, the summation in (5) over positive integer 
m > Nq takes into account the resonance (or " essentiaT ) intermediate continuum states 
corresponding to open ATI channels (to which the process of direct ATI is also possible), 
which are only supposed to give the main (predominant) contribution to HHG process 
under consideration (see also [29] for details). It is also worth noting that the time- 
independent matrix function Rn,m{^]i',(lm,v) allows for representation in the factorized 
form: 

Rn,n. (fik', Qn., V) = J ff""^ (^k', q„^, V) /if '''^ (Q™, V) dO^^ (9) 

that makes the final results transparent for interpretation. Namely, the integrand in (9) is 
a product of time-independent amplitudes of two strong-field multiphoton processes - the 
direct ATI process accompanied by absorption of m incident laser photons and emission 
of a photoelectron to intermediate continuum states of canonical momentum q^ 

/if ^^) (q^, r^) ~ 27r {Up - mu) 5„„ (^((Qm); |) ^o (q^^) (10) 

followed by subsequent process of spontaneous photorecombination and emission of a 
high-harmonic photon of frequency VL]^i and polarization ek',A' (A = 1,2): 

/f ^) (fik', q„^, V) ~ iqmpTi/^l (ej;,,,, ■ q™) 5„ (c(q™); |) F* (q„) (11) 

Finally, the high-harmonic spectrum within the proposed model is represented by a se- 
quence of equidistant (separated by a laser fundamental frequency u) discrete peaks of 
heights defined by the partial HHG amplitudes f)^ (^k') (3) and respective differential 



rates wjy v^)'- 



w 



[HHG),-, /N _ ^N 



N ^ ) ~ /c^ \3 



^(N) 



(k') = -^ D^nJ (c^) (12) 
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of emission of high-harmonic photon of discrete frequency Q^ = Nu to a fixed solid angle 
element rfOk' . 

The total ionization rates Tion {v) can be derived just from the amplitude (10) of 
direct ATI process by means of integration of its squared modulus over the angles of 



photoelectron emission as well as photoelectron energy. The integration over the final 
energy of emitted photoelectron is just reduced to the summation over the total number 
A^ of laser photons absorbed (see, also [30, 31], for details): 

r^oniv)= E ^AU,-NufldO^,B'_^{ciPN);l)\FoipN)f (13) 



N>No(v) 



In a pulsed laser field of arbitrary duration r, the total yield of photoelectrons (i.e., 
the produced photoelectron and/or ion signal that is measured in experiment) depends 
on the peak field strength Eq and laser's spatial and temporal profiles. For spatially 
homogeneous incident laser field, the amplitude E of laser field is only dependent on the 
peak field strength Eq and temporal profile function g (t) according to: 

E (t) = eE (t) sin (ut) = eE^g (t) sin (cut) (14) 

So, for example, for Gaussian shape of laser pulse intensity: 



g (t) = exp 



-2f^^)'ln2' 



(15) 



with to = r/2 as the time moment at which the laser intensity reaches a maximum during 
the laser pulse duration. Under conventional supposition that the ionization process 
is generally a considerably faster than the time period of essential change of temporal 
profile function g{t), the ionization rate (13) derived assuming time-independent laser 
field intensity can be approximately considered as still valid, but corresponding to the 
instant ionization rate. This means that the ionization rate (13) is just supposed to be 
adiabatically changing with variation of laser intensity through the laser pulse action and 
the total ionization probability for the overall time of laser pulse duration r can be found 
according to the quasistatic formula: 

Pionir) = 1 - exp I- r Tion iv {t)) dt) (16) 



where, particularly 

v{t)=E^t)/{Auj')=g'{t)E',/{Aco') (17) 

The Eq.(16) will be used further for defining the saturation ionization intensity Igat as 
the laser intensity which leads to 99% of the laser-irradiated species being ionized by the 
ending time of the laser pulse, so that Pionij) ~ 0.99. 



3. SAE-EXTENSION TO MOLECULAR HHG IN DIATOMIC MOLECULES: THE MOLEC- 
ULAR ORBITALS IN COORDINATE AND MOMENTUM SPACE. 

So far, up to this point, we didn't specify the nature of laser-irradiated system, so 
that all the relations above are supposed to be equally well applicable to arbitrary atomic 
and/or molecular species provided the SAE consideration is valid to derive the Fourier 
transform of a single-electron wavefunction of initial bound discrete state in a closed 
analytical form. Thus, according to the HHG model currently extended to molecular 
case of homonuclear diatomics, the specified properties of laser-irradiated system are 
incorporated solely through the ionization potential I^^ and Fourier transform (6) being 
the wavefunction F„ (p; Rq) of initial n-th molecular discrete state in momentum space. 
The corresponding single-electron wavefunctions, or molecular orbitals, each of a fixed 
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discrete binding energy Eq = ^Ip^^ and respective number A'^g^"^ of identical electrons, are 
the mathematical constructs used to describe the multi-electron wavefunction in molecules 
(similar to atomic orbitals in atom). According to the standard linear combination of 
atomic orbitals (LCAO) and molecular orbitals (MO) method, the wavefunction of each 
(n-th) molecular valence shell can be approximately considered as two- centered MO: 

$„ (ri; rs) = $„ (r - Ro/2; r + Ro/2) = 
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being a linear superposition of predominantly contributing (j-th) one-electron AOs 0^ (r) 
centered on each of the atomic cores and thus separated by internuclear distance Rq. Here 
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is the respective atomic orbital overlap integral. 

Because the highest-lying orbitals are responsible for chemical properties, they are of 
particular interest. For N2 molecule and the Sa^ HOMO of gerade and bonding symmetry 
commonly considered, the approximate two-centered single-electron molecular wavefunc- 
tion can be composed, for example, as a superposition of scaled hydrogen-like 2pz or Is 
atomic orbitals [31]: 
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Accordingly, for ungerade liiu inner valence shell of bonding symmetry or gerade Ivr^ 
HOMO of antibonding symmetry commonly considered in O2 molecule, the respective 
approximate two-centered single-electron molecular wavefunction can be composed, for 
example, as a superposition of scaled hydrogen-like 2px (or 2py) atomic orbitals [31]: 
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Note first that the "— " combination in the right-hand side of Eq.(18) doesn't not neces- 
sarily correspond to an antibonding valence shell (such as iTg or CTu) having a negligibly 
small electron density near the central region between the atomic nuclei in homonuclear 
diatomic. Neither the sign "+ " always corresponds to a bonding valence shell only (such 
as ag or 7r„) with a considerable electron density near the internuclear axis region. Let us 
also recall, that the "+ " superposition of contributing AOs doesn't not necessarily imply 
only gerade (viz. spatially symmetric, like ag or iTg) molecular valence shells, but it may 



be also inherent to ungerade (viz. antisymmetric, like, e.g. 7r„) bonding MOs. Accord- 
ingly, the "— " superposition may also correspond both to gerade (such as antibonding 
Tig) and ungerade (such as bonding 7r„) molecular valence shells. Thus, depending on the 
spatial parity of contributing AOs 0^" (r) , the " +" or " — " combination in the right-hand 
side of (18) is chosen to provide a required spatial and bonding symmetry of composed 
molecular valence shell. Meantime, the form of contributing AOs 0^" (r) in (20)- (23) is 
generally chosen to adequately (although, approximately) reproduce the spatial distribu- 
tion of electronic density in respective MO under consideration relative to the internuclear 
axis. In particular, an appropriate choice of predominantly contributing AOs (p^ (r) al- 
lows for approximate reproducing a considerable electron density (inherent to a valence 
shells) or negligibly small electron density (inherent to tt valence shells) nearly along the 
internuclear axis spatial region, irrespectively of bonding or antibonding symmetry. For 
example, the 3o"g HOMO in N2 is well-known as a considerably prolate along the inter- 
nuclear axis, so that respective spatial distribution of molecular electron density can be 
approximately reproduced by the linear combination (20) of two predominantly contribut- 
ing 2pz hydrogen-like atomic orbitals, which are oriented along the internuclear molecular 
axis (supposed to be coincident with spatial OZ axis). Accordingly, the Itt^ HOMO in 
O2 is well-known to have a negligible electron density along the internuclear axis and this 
feature can be adequately reproduced by a single linear combination (23) of two predom- 
inantly contributing 2p,j. (or 2py) scaled hydrogen-like atomic orbitals, which are oriented 
perpendicularly to the internuclear axis. 

We have to mention here that the " +" combination of 2pz scaled hydrogen-like atomic 
orbitals we previously used in [31] to reproduce 3o"g HOMO in N2 doesn't provide the 
even [gerade] spatial parity and bonding symmetry of the (2p) 30"^, although proved to be 
a surprisingly well working in explaining of relevant experiment [26] (which, particularly, 
shows no suppression found in strong-field ionization of iV2). Therefore, presently the "— " 
combination of 2pz states in Eq.(20) is used, which does provide the even spatial parity of 
the (2p) 3ag HOMO. However, according to the MO-SFA model [19], such a composition 
would result in destructive intramolecular interference and accordingly high suppression 
in ionization of N2 (see also [31], for details). To eliminate the mentioned deficiency, in our 
present consideration we apply a more accurate composition of the 3o"g HOMO taking into 
account some admixture of contribution from atomic s-states required mostly to provide 
a good agreement with experiment (which doesn't show a suppressed ionization in A^2)- 
Thus, the 30"^ MO is to be further approximated by the coherent superposition of few 
different MOs corresponding to separate contributions from atomic states of a specified 
orbital symmetry (viz., the scaled hydrogen- like Is, 2s and 2pz orbitals): 

<^(is2s2p)3a, (r; Ro) = ^is<^(is)3a, (r; Ro) + A2s'^(2s)3a, (r; Ro) + ^2p'^{2p)3a, (r; Ro) (24) 

with the weight coefficients Aj (j = 1,2,3) being the relative contributions {\Aj \ < 1) 
from respective atomic states and considered as variational parameters to be found from 
the equation for minimum of respective molecular binding energy, the value of which is 
put to be equal to the experimental value. Unlike this, the Ivr^ HOMO in O2 is further 
supposed to be reproduced by the single {2p) liig MO (24), which proved to be fairly 
well working [31] in interpretation of a high suppression of strong- field ionization found 
in relevant experiment [25, 26]. 

The applied velocity gauge (VG) formulation of SPA consideration of molecular 
strong-field ionization also implies analytical representation for Fourier transform F„ (q, Rq) 
corresponding to matrix element (6) of EM transition from the n-th initial molecular 



discrete state $n(r; R-o) to the final continuum state of a definite value of canonical mo- 
mentum p. For molecular valence shells Sa^ in N2 and Itt^ in O2 under consideration, the 
explicit expressions (20)- (23) allow for direct analytical calculation of respective Fourier 
transform (or molecular wavefunction in momentum space) for each of contributing one- 
electron two-centered single molecular orbitals: 

2'>kII^Pn cos (^p) sin [(pivRo) /2] 



(26) 
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F (r. Tt \ - yf^ri^r (r \ ^^^»^' ^^^ ~ ^^^ ^""^ [(P^'Ro) /'A 

nV2ip% + Klf J2[l + S2s{Ro)] 



F. . (n Ji)- xlJ^^Ci, ^ 2^/^/.^/^ cos[(p^-Ro)/2] 
F(.)3..(p.,Ro)-VA.e ^^^^ \^^.^^,.Y J2[I^SUR.)] ^''^ 



2' kI/'pn sin (Op) sin [(p^-Ro) /2] 



w.(P..R.)--vA'r-'c,.„, ^^^^^^, /^^T5;:(7w 



with the polar angle 6p of photoelectron emission relative to internuclear molecular axis 
(cos(6'p) = (p ■ Ro) /pRo) and A/'f'"«) = 2 and N^^"":^^ = 4. Here k„ = z]"V«i , whereas 
Zj is the effective charge corresponding to " ejfective" long-range Coulomb model bind- 
ing potential of respective residual molecular ion, while Oj = ja^ is j-th Bohr orbital 
radius of respective atomic orbital contributing to n-th initial molecular discrete state of 

molecular binding energy €q = — /t^/2 = — J^"-' = — (Zj /aj) /2 under consideration 

(viz., e'i'''^ {N2) ^ -15.58 eV, ej,^"^^ (O2) ^ -12.13 eV). The analytical expressions for 
respective atomic orbital overlap integrals (19) are presented in [31], moreover, due to 

the VG version of SFA applied, the correction factor C (/t„) = f2/t„/^"Y_Ej " is also 
introduced to matrix elements (25)-(28) to incorporate the long-range electron-molecular 
ion Coulomb interaction in the final continuum Volkov states into account [19]. 

Let's note first that the angular dependence of the (2p) ScTg molecular wavefunction 
in momentum space (25) clearly indicates that a separate contribution from 2pz states 
alone is highly resistant to ionization along the direction perpendicular to the internu- 
clear molecular axis (cos (^p) ~ 0). The respective angular dependence on the angle 6p of 
photoelectron emission relative to the internuclaer axis is well seen in Fig. la to be a consid- 
erably prolate along the internuclear axis. Thus, the photoelectron emission from (2p) Sag 
MO is profoundly dominated along the internuclear molecular axis (cos (^p) ~ 1). similar 
to observed earlier in [31] for incorrect (the "+") combination of two 2pz states. This, 
particularly, illustrates that the spatial (gerade or ungerade) symmetry of the (2p) Sag 
MO doesn't essentially affect the angular dependence (viz., nodal plane along the molec- 
ular axis) of respective Fourier transform. Owing to the presence of cos (^p) factor (see 
Eq.(25)), both the "— " (currently used) and the "+" (earlier used in [31]) combination 
of 2pz atomic states suggest a predominant photoelectron emission along the internuclear 
axis (see also Fig. la presented in [31], for comparison). The sign of the combination 
mostly alters the absolute value of local maxima in respective Fourier transform due to 
different {sine or cosine) trigonometric factor arising due to intramolecular interference 
of contributions of ionization from two atomic centers [19] separated by the internu- 
clear distance Rq. Since the generalized Bessel functions in Eqs.(10)-(11) are commonly 
maximized when photoelectron momentum p is parallel to the laser polarization e, the 
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ionization from {2p) Sag HOMO composed as a single combination of 2pz states alone 
is thus expected to be predominant when A'^2 molecular axis is aligned along the laser 
field polarization (Ro||e). The latter conclusion (also previously made in [31] for {2p) 3crg 
HOMO of ungerade symmetry) seems to be natural and consistent with relevant experi- 
ment [32] and earlier alternative VG-SFA calculation [19] (see Fig. 7b presented therein) 
using pure numerical self-consistent Hartree-Fock (HF) procedure of 3o"gHOMO composi- 
tion taking the contribution from atomic s-states into account. However, this statement 
is in a contradiction to respective VG-SFA results of calculation [33] applying the same 
VG-SFA approach and numerical HF-based procedure of Sag composition as proposed 
earlier in Ref. [19]. Contrary to [19], the strong-field ionization of N2 was found in [33] 
to be predominant for perpendicular orientation of internuclear axis relative to the laser 
field polarization. 

To clarify the reason of the contradiction. Figures lb- Id demonstrate the angular 
dependence of molecular wavefunction in momentum space calculated for possible differ- 
ent 3ag compositions (5), either with or without taking a comparable contribution from 
atomic s-states into account. The relative (weight) coefficients Aj at Is , 2s and 2pz states 
under consideration proved to have comparable values for Sag composition presented in 
Fig.l, although not necessarily of the same signs. For example, Ais ~ — ^2p ~ 1/a/2 for 
{ls2p)Sag, whereas Au ~ — A2S — — ^2p ~ 1/a/3 for (Is2s2p) 3crc, composition, so that 
the weight coefficients at Is and 2pz states (as well as Is and 2s states) are found to have 
always opposite signs. Moreover, only the combinations with Ais and A2p (as well as Ais 
and A2s) of opposite signs proved to provide no suppression in N2 ionization revealed 
in experiment [26]. To summarize, the presented Fourier transforms corresponding to 
{ls2p)Sag and/or {ls2s2p)Sag valence shells clearly suggest a predominant photoelec- 
tron emission from Sag HOMO along the internuclear axis, contrary to suggested in [33]. 
Our Fig.l also demonstrates that such difference may only arise owing to relative contri- 
bution from s-states, which seems to be somehow overestimated (compared to p-states) 
within the particular (12s7p)/[6s4p] basis chosen in [33] for numerical HF-based proce- 
dure of Sag composition in coordinate space (see also [34] for details). In the meantime, 
a composition overestimating the relative contribution from s-states doesn't seem to be 
suitable to Sag molecular state, which is known to have a more complex structure surely 
dominated hj p-states in coordinate space [17]. This also becomes especially evident ow- 
ing to Fig. 12 in Ref. [35] displaying the correct orbital shape for the Sag wavefunction in 
A^2j which obviously contrasts with the resulting Sag state of substantially different shape 
and detailed structure suggested in [33] (cf., for example, respective Fig.lc presented in 
the last reference therein). 

For laser-irradiated atomic counterparts of diatomics under consideration - such as 
Ar (for A^2 and/or H2) and Xe (for O2) with Sp and 5p outermost atomic valence shells, 
respectively - the initial atomic ground state can be also approximately reproduced by 
corresponding one-electron hydrogen-like atomic orbital. The analytical expression for 
respective atomic Fourier transforms F^p^'" ^^ (p) of Spz, Sp^ and Spy hydrogen-like AOs 
(e.g., in Ar) corresponding to different possible values of associated magnetic quantum 
number m {0 - for Spz and ±1 - for Sp^ ± iSpy) has the form: 



-.{3p^±i3py) , ^ _ iC (/€„) 2'^^/2 k}^/'^pn ( ± sin {d) exp {±iif) 



TTA 



Here the angles 6 and ip are the polar and azimuthal angles of photoelectron emission 
with respect to the incident field polarization (cos {6) = (e ■ p) /p), so that the angular 
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factor in the right-hand side of (29) corresponds to the case of either Spx^iSpy (the upper 
part), or 3pz (the lower part) AO. The "effective" charge Z^Vr of pure Coulomb model 
binding potential corresponds to the experimental value of atomic ionization potential 

Jp = kI/2 = (Z^^l/ag) ^2. {Ip ~ 15.75 eV - for Ar and Ip ^ 12.13 eV - for Xe). 

The proposed strong-field model is also equally well suitable to describe molecular 
HHG process in arbitrary molecular valence shell including inner ones closest to the 
HOMO (although, of higher binding energy, different orbital symmetry and number of 
identical electrons). The contributions from different (outermost and inner) valence shells 
are to be added coherently in molecular HHG, so that total amplitude of emission of A^-th 
harmonic is just the sum of respective partial amplitudes corresponding to contribution 
from each molecular valence shell under consideration. In particular, the overall A^-th 
partial amplitude (4) of harmonic emission in diatomic N2 can be represented as consisting 
of contributions from three highest molecular valence shells - Sag (J^^'^"-' ~ 15.58 eV), Itt^ 
(^jiiTTu) ^ iQQQ g-j/) and 2a„ (J^^a^) ^ ^g 73 ^yy^ 

/f^^ (fikO = fS'''^ (aO + A''^"^ (fikO + A''^"^ (aO (30) 

whereas, for HHG in diatomic O2, accordingly, as contribution from three different highest 
molecular valence shells - lUg {ij^^^^^ ^ 12.07 eV), l7r„ (J^^^-) ^ 16.26 eV) and Sag 

/^^ (ao = fi^^'^ (ao + A'""^ (ao + fS'''^ (ao (31) 

Under the conditions of equivalent atomic centers and the atomic or molecular specie 
occupies the same state before and after the passage of the laser pulse, only odd harmonics 
of the fundamental laser frequency of linear polarization are emitted [29]. The harmonic 
emission rates are to be currently calculated only for high-order harmonics of polarization 
ek',A', for which the emission rate becomes maximal, that is parallel to polarization e of 
incident laser field. The model-related background analytical expressions (3)-(12) and 
(25)-(28) do imply fully SD consideration of the process (particularly, they are available 
for arbitrary orientation of internuclear molecular axis with respect to the laser polar- 
ization e), nonetheless, the molecular axis is supposed further to be strongly aligned 
along the polarization of incident laser field (i.e. (e ■ Rq) ~ cos (6) = ) throughout 
our present consideration, unless stated otherwise. The associated orientation effects in 
molecular harmonic generation are thus beyond the scope of our present consideration, 
which is mostly focused on behavior of molecular harmonic spectra depending on equally 
important problem parameters, such as internuclear separation Rq and properties of a 
specified laser-irradiated diatomic (such as orbital symmetry, composition and structure 
of molecular valence shells) as well as incident laser parameters (such as laser intensity, 
pulse duration and temporal profile). 



4. CALCULATION, NUMERICAL RESULTS AND DISCUSSION. 

When studying the strong-field phenomena, the major difficulty with molecules is 
to deal with both the electronic and the nuclear degrees of freedom (i.e., related to the 
molecular vibrational and rotational motion). Nonetheless, the including of the nuclei 
motion into a proper consideration is generally feasible within frame of the currently ap- 
plied approach, but this additionally requires of preliminary knowing available reliable 
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data (previously calculated or measured in experiment) related to dependence of molec- 
ular binding energy on internuclear separation Rq. To the best of our knowledge, the 
latter dependence is reliably known only for the simplest diatomics {H2 or H2), solely for 
which the Rq -dependent behavior of molecular HHG can be numerically studied within 
the currently proposed model. For all other diatomics under consideration (A^2 and O2), 
the nuclei motion will be further completely ignored and the Born-Oppenheimer approx- 
imation is supposed to be valid. For this reason we start our numerical study from H2 
molecular ion, which is sufficiently well explored in previous theoretical and experimental 
studies (see [36], for a review). Beginning from the first HHG calculations [37, 38], this 
molecular specie, owing to its structural simplicity, became the object of a large number 
of other theoretical and experimental investigations (e.g., [8, 10-12]), the results of which 
can be used for direct comparison. 

Fig. 2 displays our results for high-harmonic spectra of H2 calculated as harmonic 
emission rates (13) compared to respective atomic counterparts of nearly identical binding 
energy under conditions considered earlier in [10]. Like in [10], the results presented in 
Fig. 2a were calculated assuming /?o-iiidependent fixed molecular binding energy Eq (viz., 
I^ol = 0.4 a.u. for any internuclear separation Rq), which is supposed to be identical to 
that of model atomic "counterpart" bound to model zero-range potential (ZRP). Since 
the molecular binding energy is known to be i?o-dependent, Fig.2b also presents the 
analogous high-harmonic spectra of H2 calculated under the same laser pulse and Rq- 
dependent molecular binding energy Eq (Rq) corresponding to initial lag ground state. 
Also, unlike Fig. 2a, the high-harmonic spectrum of respective atomic "counterpart" was 
calculated assuming the hydrogen-like IS* discrete state of atomic binding energy eq ~ 
— 15.86 eV nearly identical to the molecular binding energy Eq ~ —0.583 a.u. of H2 (lo'g) 
corresponding to Rq = 12 a.u. The molecular harmonic spectra presented in Fig. 2 show 
the clear signs of the so-called " low-energy hump" similar to that revealed in [10] within a 
low-frequency domain, which is enormously enhanced for two-centered case as compared 
to the one-centered ("atomic") case also presented in Fig. 2. The origin of this "hump" 
arising in molecular high-harmonic spectra was well explained earlier [10] in terms of the 
rescattering mechanism and, particularly, by the existence of two different pathways for 
the initially released electron to recombine: with the center it started from as well as 
with the other. Namely, the "hump" corresponds to the case when the initially released 
electron moves directly from one center to the other within a time interval that is short 
compared with the period of the field, moreover, this happens around the time where the 
field is near its maximum. Under these conditions the energy gain, which the released 
electron acquires from the driving laser field when moving along the shortest way from 
one atomic center to other, can be approximately evaluated just by Eei ~ ERq that 
corresponds the position of the top harmonic N^ ~ ERq/ 00 of the "hump". According to 
this quasiclassical interpretation, the position of its top harmonic should be noticeably 
shifted toward a higher frequency part for a larger internuclear separation Rq. This 
particular behavior is fairly well reproduced in Fig. 3 by other high-harmonic spectra of 
H2 calculated for a different incident laser pulse and various larger values of Rq, for which 
Fig. 3 also shows the signs of "split" arising within the "low-energy hump". 

We also note the position of high-frequency plateau cutoff calculated for relatively 
large internuclear separation and presented in Figures 2 and 3 noticeably deviates from 
what is prescribed by the well known rule " Ip + 3.17f/p" for standard harmonic cutoff 
frequency. So, according to this rule, one would expect for the position Nc of the cutoff 
frequency to be nearly at the harmonic order Nc ~ 43 in Fig.2a or Nc ~ 49 in Fig.2b, that 
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is quite consistent with those figures suggest for atomic case. Meanwhile, the respective 
molecular harmonic spectra calculated for Rq = 12 a.u. clearly demonstrate a noticeably 
different position for the cutoff frequency, viz., at harmonic order Nc ~ 51 in Fig.2a or 
Nq ~ 55 in Fig. 2b. Moreover, there is also a pronounced and well-defined " cutoff splif 
visible within the plateau cutoff domain of molecular harmonic spectra presented in Figs. 2, 
3. This "cutoff split" is similar to the split revealed earlier [10] in the cutoff domain of 
H2 high-harmonic spectra calculated for sufficiently large internuclear separation Rq. 
The origin of the "cutoff split" can be also well explained in terms of the rescattering 
mechanism and, particularly, attributed to the interference of two other (longest, but 
topologically different) pathways for the released electron moving from the one atomic 
center to the other to recombine with (see, also [10], for details). The respective separation 
in energy the released electron gains from the driving laser field is about at the position 
of high-frequency plateau cutoff (viz., nearly around the mentioned above harmonic order 
A^ = 47 in Fig. 2a and/or iV = 51 in Fig. 2b), so that the two maxima alongside completely 
mask the actual cutoff position as the distance Rq increases. 

As noted above, unlike the atomic case, the cutoff position Nq in calculated molecular 
harmonic spectra proved to be not fixed yet even for the fixed parameters of incident laser 
pulse, but also substantially dependent on internuclear separation. Such i?o-dependent 
behavior of Nc (Ro) is well illustrated by Fig. 3, which represents H2 harmonic spectra 
calculated for different fixed values of Ro under the same laser pulse. In particular. Fig. 3 
demonstrates the deviation of the harmonic cutoff position from the well known rule 
"/p -|- S.nUp" , even if the corresponding decrease of molecular ionization potential Ip for 
larger Rq is taken into complete account. So, in Fig.3a the position of high-frequency 
cutoff for relatively small Rq = 4 a.u. is in a full accordance with the aforementioned 
rule for standard cutoff corresponding to the harmonic order Nc ~ 73. This value is 
accordingly and noticeably smaller compared to A'",^ ^ 80 prescribed for larger ioniza- 
tion potential Ip = 30.0 eV corresponding to the equilibrium internuclear separation 
Rq = 2 a.u. Thus, the calculated position of harmonic cutoff initially moves toward a 
lower frequency domain of spectrum due to corresponding decrease of molecular ioniza- 
tion potential Ip as Rq increases. However, under further expansion and corresponding 
decrease of molecular ionization potential, the cutoff frequency is, nevertheless, getting to 
move back toward a higher frequency domain, contrary to the "standard cutoff" rule. So, 
for example, in molecular harmonic spectra calculated for larger internuclear separation 
and presented in Fig.3a, one can identify the cutoff position located at harmonic order 
Nc {Rq = 15 a.u.) ~ 79 and Nc {Rq = 30 a.u.) ~ 89 respectively. These two values are 
noticeably larger than Nc ~ 70 prescribed by the mentioned " standard cutoff" rule based 
on correspondingly decreased molecular ionization potential in H2 taken into account. 
The physical mechanism underlying this quite an opposite tendency of increasing the har- 
monic cutoff frequency for sufficiently large Rq increased can be again well understood in 
terms of semiclassical "recollision" picture invoked above for interpretation of the "cutoff 
split". Namely, when Rq increases, the released electron acquires an additional energy 
gain from the driving laser field when moving along the longer ways from one atomic cen- 
ter to other. Due to the mentioned additional energy gain followed by recombination the 
acquired electron energy is released as a photon of a considerably higher frequency that 
gives rise to increasing of the harmonic cutoff frequency. Beginning from sufficiently large 
internuclear separation (over about Rq = 10 a.u.), this mechanism becomes predominant, 
insomuch that the related additional energy gain can compensate or even overcome the 
competing decrease of standard cutoff frequency due to diminished molecular ionization 
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potential Ip for larger Rq. Such a complicated -Rq- dependent behavior of harmonic cutoff 
position revealed in molecular high-harmonic spectra under discussion seems to be very 
promising for various applications. To the best of our knowledge, the dissociating di- 
atomics with large internuclear separations were predicted first in [39] to be capable to 
produce high harmonics of frequency far beyond the "standard" (atomic) cutoff, however, 
this idea has not been confirmed yet by experiments. Particularly, when the nuclei are 
separated by a distance equal to odd multiples of vrao (where a^ = {E /u"^) is the classical 
amplitude of electron's excursion in laser field), the cutoff limit was shown to extend 
beyond Ip + 3.17[/p up to Ip + 8f/p since the limit of 8f/p is the maximum kinetic energy 
reachable by the electron during its classical excursion from one atomic center to other 
(see also [10]). The results presented in Fig. 3 generally confirm such a prediction and 
indicate unambiguously that, quite a noticeable extension of high-frequency plateau cut- 
off is to be present in molecular harmonic spectra for very large internuclear separation. 
Although, due to relatively small values of Rq used under calculations compared to tiuq 
^110 a.M., the calculated extension of high-frequency plateau is not very impressive. 

On the other hand, the values of Rq considered in Figures 2, 3 proved to be sufficiently 
large to confirm the result of previous research [10, 38] showing that molecular harmonic 
spectra do not obey the "Jp -|- 3.17f/p" rule yet for the position of high-frequency cutoff. 
Moreover, unlike atomic high-harmonic spectra, the molecular harmonic cutoff is split 
due to intramolecular interference and this "cutoff split" well seen to be a noticeably 
broader and more pronounced for larger Rq. Also, like in [10], the currently calculated 
H2 harmonic spectra in Fig. 3 clearly demonstrate a significant enhancement (up to several 
orders of magnitude) in harmonic emission rate as internuclear separation Rq increases. 
This is not surprising for diatomics since the respective molecular binding energy \eQ {Rq)\ 
is to be smaller for larger Rq. The similar enhancement of harmonic emission due to 
decreased ionization potential was initially revealed in atomic high-harmonic spectra (see, 
e.g., Fig.5 and Fig. 10 in [40] as well as references to relevant experiments cited therein). 
This also allows for quite a clear explanation, namely, the laser-irradiated species having 
a smaller ionization potential are generally more polarizable, so that their response to 
the incident laser field is expected to be a noticeably stronger than in harder ionized 
species. Our results calculated for II2 and presented in Fig. 3 also confirm that the similar 
underlying physical mechanism is responsible for the enhancement of harmonic emission 
rate in diatomics for larger internuclear separation as well as the transient enhancement 
of the high-harmonic conversion efficiency recently found as exceeding the atomic one 
in expanding molecules [12, 14]. It is also very interesting that Fig.3 still demonstrates 
noticeable differences in details of two high-harmonic spectra of i/^ corresponding to 
different lag and la^ initial molecular states and calculated for large Rq = 30 a.u., at 
which generic structural features (e.g., such as orbital symmetry) of particular laser- 
irradiated diatomic are expected to have a very minor importance. 

Another remarkable feature well recognized in molecular harmonic spectra presented 
in Fig. 3a is a sufficiently pronounced local minimum clearly visible within the high- 
frequency plateau region and located respectively at A^ = 55 for Rq = 4 a.u., N = 53 for 
Rq = 15 a.u. and A^ = 49 for Rq = 30 a.u. The minima seem to be similar to those revealed 
earlier in numerically calculated molecular high-harmonic spectra [11] and predicted to 
arise in harmonic spectra of diatomics with HOMO of bonding symmetry due to destruc- 
tive intramolecular interference of photoelectron emission from separated atomic centers 
to intermediate continuum states. Like in [11], the interference- related minimum in Fig. 3a 
shows a regular dependence of its position on internuclear separation and, particularly. 
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moves toward a lower frequency domain with increase of Rq. Such i?o-dependent behavior 
is generally in accordance with the semiclassical (although, quite an ad hoc) interpretation 
given in [11] and related rule N^^^ = Ip/u + tt^/ {2ujRq) ~ tt^/ {2ujRI) suggested for 
location of the minimum. However, the harmonic spectra calculated under the conditions 
of [11] and presented in Fig. 3a doesn't seem to confirm the above prediction for location 
of the minimum, but rather suggest a different (although, quite a raw) approximate rule 
-^min (-^p + Up)/uj + 7r^/ {2ujR^ . The presently suggested empirical rule seems to 

be similar to revealed earlier for the interference-related minimum arising in calculated 
molecular photoelectron spectra [10, 31] due to direct ATI in diatomics with HOMO of 
bonding symmetry. Therefore, the interference-related phenomena arising in molecular 
harmonic spectra should also allow for transparent interpretation similar to that given 
for molecular ATI [31]. In our opinion, this also remarkably demonstrates a close and 
intimate interrelation of molecular ATI and HHG processes. 

Since the harmonic generation process is strongly related to ionization, one should 
also expect the enhanced harmonic emission rates in laser-irradiated species for which 
the ionization rate is also enhanced for some reason. On the other hand, for laser pulses 
of sufficiently high peak intensity / and long pulse duration r, the enhanced ionization 
leads to a faster saturation of ionization and exhaustion of initial state of laser-irradiated 
system, which, after reaching the saturation of ionization. is unable to efficiently produce 
high-order harmonics anymore. This gives rise to quite a natural restriction imposed 
on standard cutoff frequency Lp + 3.17f/p (with ponderomotive energy Up ~ /), which 
thus cannot be extended yet under further increase of incident laser intensity / over 
the saturation intensity I sat- Hence, the laser- irradiated species, such as H2 and O2; 
showing a suppression in ionization rate (compared to counterpart systems of nearly 
identical ionization potential Ip) are expected to produce harmonics of considerably higher 
frequency corresponding to a longer high-frequency plateau with an extended cutoff due 
to a noticeably higher laser intensity Isat that saturates ionization. The extension in high- 
harmonic spectra of O2 relative to atomic Xe under respective laser saturation intensity 
was recently observed in experiment [13] and also partly reproduced in related MO-ADK 
calculations [13] based on the Lewenstein atomic HHG model [22]. 

In order to assess a possible extension of high-frequency plateau in molecular har- 
monic spectra of more complex diatomics, such as N2 and O2, compared to respective 
atomic counterparts, one has to evaluate first the respective saturation laser intensity for 
each of laser-irradiated species under consideration. Fig.4 displays our results for total 
probability of ionization of O2, Xe, A^2 and Ar calculated under conditions of experiment 
[13] according to formula (16) extended to the diatomic case based on earlier developed 
SFA-LCAO model [31] as described in Sec. 3 above. The respective saturation laser in- 
tensity Isat is well identified in Fig.4. and defined as the intensity at which about 99% of 
laser-irradiated species being ionized by the end of the laser pulse duration. The results 
have been summarized in the Table I, which, for the convenience of direct comparison, 
also represents the respective results of [13] calculated according to MO-ADK theory 

Table I 
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Specie 


N2 


Ar 


O2 


Xe 


HOMO 
AO 


(ls2s2j9)3CTg 


3p6 


(2p) iTT, 


5p6 


i?o(A) 


1.1 




1.21 




/p {eV) 


15.58 


15.75 


12.07 


12.13 


hat (W/cm') 

present 
calculation 


3.75 ■ 10^^ 


4.4 ■ 10^^ 


2.5 ■ 10^^ 


1.7 -10^^ 


Isat {W/cm') 

MO - ADK 

calculation[13] 


3.71 ■ 10^^ 


4.4 • 10^^ 


3.01 ■ 10^^ 


1.75-10^^ 



The table demonstrates a fairly good consistence of our currently calculated results with 
the prediction of MO- ADK theory, except the saturation intensity calculated for O2. 
The deviation is explained by a different value of molecular ionization potential, viz. 
-^p (O2) = 12.56 eV used in [13] (see, e.g.. Table I presented therein), which is a slightly 
higher compared to the currently used experimental value Jp (O2) = 12.07 eV [26]. This 
inevitably leads to an overestimated high suppression in ionization of O2 relative to Xe, 
thereby giving rise to a noticeably overestimated value of the saturation laser intensity 
3.01 ■ 10^'^W/cm'^ as compared to 2.5 ■ lO^^W/cm'^ currently calculated. 

The high-harmonic spectra of O2 presented in Fig. 5 demonstrate a noticeable exten- 
sion compared to the spectrum of Xe calculated at respectively lower laser saturation 
intensity previously found (Fig. 5b). At the same time, the extension of molecular high- 
frequency plateau and respective cutoff frequency currently calculated for O2 (relative to 
Xe) noticeably deviates from the similar extension found in [13] for different saturation 
intensities calculated using the MO-ADK model. Particularly, the harmonic order Nc 
for cutoff position in presently calculated harmonic spectrum of O2 is about 39 (versus 
calculated 43 and measured 53, both reported in [13]) provided the contribution from the 
liTg HOMO is only taken into account, like in [13]. However, the total harmonic spectrum 
of O2 presented in Fig. 5a and corresponding to the coherent contribution (31) from the 
(2p) iTTg HOMO and two closest {2p) Itt^ and {Is2s2p) 3(jg inner shells additionally taken 
into account demonstrate a noticeably longer high-frequency plateau of cutoff located 
at Nq (O2) = 45 and thus extended relative to Xe, for which the cutoff is located at 
Nc {Xe) = 29. Note also that the additional extension clearly visible in Fig. 5a is mostly 
due to the contribution from the inner Sag valence shell, whose separate contribution 
within the cutoff domain is comparable with that from the (2p) l-Kg HOMO and always 
prevailing with respect to that from the {2p) ItTu, despite the inner Sag has the higher 
ionization potential than liTg or Itt^. This noted feature seems to be quite a natural 
for strongly aligned diatomics under consideration and can be explained by an enhanced 
ionization (initial release of electron to intermediate continuum states) from the 3a g MO, 
which is a considerably prolate along the internuclear axis. Unlike a orbitals, vr molecu- 
lar orbitals (such as Itt^ HOMO and inner l7r„ in O2) are predominantly oriented along 
the direction perpendicular to the internuclear axis, so that the photoelectron emission 
(initial release of electron to intermediate continuum states) from tt molecular orbitals in 
aligned diatomics is a considerably suppressed along the internuclear axis that leads to a 
relatively smaller harmonic emission rate (see also the relevant discussion in [31]). 

The high-harmonic spectra of A^2 and its atomic counterpart Ar calculated for respec- 
tive saturation intensities currently found and shown in Table I, are presented in Fig.6. 
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In particular, Fig. 6b demonstrates some minor extension of cutoff in harmonic spectrum 
of Ar relative to N2 due to a slightly enhanced ionization rate of N2 observed in exper- 
iments [25, 26] and related lower laser saturation intensity compared to Ar. The cutoff 
position in currently calculated harmonic spectrum of N2 corresponding to the separate 
contribution from the {Is2s2p) ScTg HOMO (only taken into account in [13]) is located at 
harmonic order about 55 that is well consistent with that calculated in [13]. Meanwhile, 
the total harmonic spectrum of N2 presented in Fig. 6a and calculated taking the coherent 
contribution (30)from two closest inner {2p) Ivr^ and (2s) 2o"„ valence shells into account 
demonstrates a bit longer plateau of cutoff position at harmonic order Nq {N2) = 61 
that is in a reasonably good accordance with experimental value measured in [13]. This 
noted extension of high-frequency plateau cutoff is also clearly visible to arise mostly due 
to 2au valence shell, whose separate contribution within the cutoff region is comparable 
with that from the Sag HOMO and always prevailing over that from Ivr^, even though 
the ionization potential of 20"^ is higher. Thus, the additional extension due to the con- 
tribution from inner valence shells is quite similar to that we observed earlier in the total 
high-harmonic spectrum of O2 presented in Fig. 5a. This probably might also explain the 
difference between the positions found for harmonic cutoff in high-harmonic spectra of 
N2 calculated and observed in [13]. 

All the results presented in Fig. 5 and Fig. 6 have been summarized in Table H below, 
for the convenience of direct comparison with experiment [13]. 

Table II 



Specie (AO, MO) 


/p (eV) 


Nc 
presently 
calculated 


Nc 

MO - ADK based 

calculation [13] 


Nc 

measured in 

experiment [13] 


N2[ils2s2p)3ag] 


15.58 


55 


55 




N2 [{2p) IttJ 


16.96 


57 






N2 [(2s) 2a J 


18.73 


57 






N2 {3cxg + l7r„ + 2(Tu) 




61 




57 


Ar (3/) 


15.75 


67 


65 


63 












O2 {2P) lUg] 


12.07 


39 


45 




O2 (2p)l7r„ 


16.26 


43 






O2 {Is2s2p)3ag 


18.18 


43 






O2 (iTTg + l7r„ + 3o-g) 




45 




53 


Xe {5p^) 


12.13 


29 


29 


29 



As a final remark to the results presented in Table H, we have to comment the respective 
HHG results for N2 and Ar reported in [13] and based on MO-ADK calculation of laser 
saturation intensity (see Table I in Ref. [13]). Namely,the latter results seem to be at least 
confusing and contradicting to the MO-ADK results presented ibidem in Fig. 4, which un- 
ambiguously shows a considerable suppression for A^2^ ion signal compared to Ar^ for the 
shorter laser pulse of duration r = 30 fs. Under these conditions close to experiment [13] 
corresponding to r = 25 fs, the MO-ADK model suggests the calculated ratio N2 /Ar^ 
for ion signals from A^2*^ relative to Ar^ to be nearly around 0.5 (see also Fig.8 in [23]) that 
is consistent with relevant experiment [25], at least within the strong-field domain where 
the tunneling ionization theory is only valid. This, however, implies that the MO-ADK 
based calculation would result in an accordingly higher saturation intensity rather for N2 
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than for Ar, similarly to the suppressed ionization in O2 relative to Xe. Contrary to such 
a natural conclusion the calculated laser saturation intensity for N2 was identified in [13] 
to be a noticeably lower (viz., 3.71 ■ 10""^^ W/cm'^ versus 4.4 ■ 10^^ W/cm'^ reported for 
Ar) that seems to be self-contradicting and thus makes a fairly good agreement of the 
calculated and measured harmonic cutoff frequency of respective high-harmonic spectra 
observed in [13] quite a questionable and controversial. 



5. CONCLUSION. 

The strong-field process of high-order harmonic generation was considered and stud- 
ied numerically in a number of laser- irradiated homo nuclear diatomics (namely, H2 , 
N2 and O2) compared to their atomic counterparts of nearly identical ionization poten- 
tial. The proposed strong-field molecular HHG model is essentially based on quantum- 
mechanical SFA-based approach developed earlier and currently extended to molecular 
case using the standard LCAO-MO method applied for adequate reproducing the wave- 
function of initial molecular discrete state. The latter wavefunction is approximately 
considered as a single- electron and two-centered molecular orbital (MO) consisting of 
coherent superposition of few predominantly contributing scaled hydrogen-like atomic 
orbitals (AO) centered on each of the atomic cores and thus separated by internuclear 
distance Rq. The form of contributing AOs is generally different for various molecular va- 
lence shells (such as ag , iiu and vr^) of different orbital and bonding symmetry. Wherever 
appropriate, the comparable contributions from other (inner) molecular valence shells of 
a larger binding energy and different orbital symmetry are also incorporated. This, par- 
ticularly, allows for fully 3D consideration of molecular HHG process as well as deriving 
the final expressions for amplitude and rate of harmonic emission rate in closed analytical 
form available for further numerical calculations. 

The calculated molecular high-harmonic spectra were found to exhibit a general shape 
similar to atomic ones, with extended high-frequency plateau ending by a well defined 
harmonic cutoff. Although, the detailed structure, such as the extent of high-frequency 
plateau and respective harmonic emission rates demonstrate a number of noticeable and 
remarkable differences from respective atomic harmonic spectra. These distinct differ- 
ences, such as the "low-frequency hump" and " cutojf splif arising due to two-centered 
nature of molecular bound state, are ascertained to be strongly dependent on the internu- 
clear separation Rq. For sufficiently large Rq the molecular harmonic spectra demonstrate 
a plateau of a noticeably longer extent and significant enhancement in harmonic emission 
rates compared to spectra corresponding to equilibrium internuclear separation and cal- 
culated under the same laser pulses. In addition, our SEA calculations seem to confirm 
the existence of a broad and pronounced minimum predicted to arise in harmonic spectra 
of diatomics with HOMO of bonding symmetry due to destructive intramolecular interfer- 
ence of contribution from photoelectrons of large momenta emitted from separate atomic 
centers [11] to intermediate continuum states. Namely, our results presented in Fig. 3a 
seem to reproduce i?o-dependent behavior of the interference-related minimum found in 
ID and/or 2D TDSE calculations [11], although also suggest a bit different (compared 
to [11]) approximate rule for location of this minimum in molecular harmonic spectra. 

Finally, the average harmonic emission rate proved to be very sensitive to orbital 
symmetry of molecular orbital and its predominant orientation with respect to the in- 
ternuclear molecular axis, insomuch that for some group of harmonics the contribution 
from inner molecular shell of higher binding energy (but, different orbital symmetry) can 
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be comparable (or even predominant within the plateau cutoff region) relative to that 
from the HOMO. Also, due to a high suppression revealed in ionization rate for O2 rela- 
tive to Xe atom, the harmonic spectrum calculated for O2 was found to demonstrate a 
high-frequency plateau of a considerably longer extent (and accordingly higher harmonic 
cutoff frequency) than in spectrum of Xe atom, for which the saturation laser intensity 
was calculated to be a noticeably smaller. To conclude, the numerically studied behavior 
of harmonic spectra calculated for O2, N2 and H2 depending on the problem parameters 
proved to be in a reasonable or fairly good accordance with relevant data calculated by 
earlier authors and observed in molecular HHG experiments. 
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FIGURE CAPTIONS 

1. Fig.l. (Color online). The squared module Fs^-g (pat, Rq) of molecular wavefunction in 
momentum space (Fourier transform) for N2 vs the angle d-p (relative to the internuclear 
molecular axis lined up in the horizontal direction) of photoelectron emission from the 
3(Tg HOMO of different composition: (a) (2p) 3o"g MO composed from 2pz states alone, 
(b) {ls)?)ag MO composed from Is states alone, (c) (ls2p) Sa^ MO composed from Is 
and 2pz states, (d) (Is2s2p) So-g MO composed from Is, 2s and 2pz states. Those angular 
dependencies are presented for photoelectrons only emitted along the laser field polariza- 
tion due to absorption of different number A^ of laser photons beginning from minimum 
one A'^o = 18 required for ionization of the 3(Tg HOMO by Ti:sapphire incident laser field 
of the wavelength A = 800 nm {huj = 1.55 eV) and fixed intensity / = 2 ■ lO^^W/cm'^ . 

2. Fig. 2. (Color online). The differential harmonic emission rates w}^ (k') [Eq.(13)] in 
strongly aligned H2 vs the harmonic order A^ calculated under conditions of Ref. [10] 
for incident laser field of the wavelength A ~ 1064 nm {hu = 1.172 eV) and intensity 
/ = 1.19 ■ 10^^ W/cm'^ (so that r] ~ 10.58). The presented molecular high-harmonic 
spectra (filled circles) are for the internuclear separation Rq = 12 a.u. assuming (a) 
the fixed (independent on Rq) binding energy Eq (|eo| = 0.4 a.u.) (b) i?o-dependent 
molecular binding energy Eq (Rq) corresponding to (Is) IcXg valence shell. The respective 
spectra of model counterpart "atom" (open circles) of the same ionization potential Ip 
were calculated assuming (a) the model zero-range binding potential (ZRP), (b) pure 
Coulomb binding potential. The results of Fig.2a are to be compared with relevant ones 
calculated within a different HHG model [10] (see Fig.4 presented therein). 

3. Fig. 3. (Color online). The differential harmonic emission rates w^ (k') [Eq.(13)] vs 
the harmonic order A^ calculated for various internuclear separations Rq in H2 strongly 
aligned along the polarization of Tiisapphire laser field of the wavelength A ~ 800 nm 
{huj = 1.55 eV) and intensity / = 5 ■ 10^'^W/cm'^. The calculated harmonic spectra are 
presented for two different bonding symmetry of initial molecular state: (a) bonding lag 
molecular state, (b) antibonding la^ molecular state. The molecular binding energy eq is 
supposed to be dependent on internuclear separation _Ro, so that the dependence Eq (Rq) 
is to be different for la^ and lag molecular orbitals. Beginning from Rq = 10 a.u. the 
dependence Eq (Rq) becomes almost the same for lcr„ and lag valence shells and, to a fairly 
good accuracy, can be approximately reproduced by the equation Eq (Rq) ~ —0.5 — 1/Rq. 
These results are to be compared with relevant ones presented in Ref. [11] (see, e.g., Fig.4 
therein). 

4. Fig.4. (Color online). The temporal dependence of total ionization probability Pioni'^) 
[Eq.(17)] calculated under conditions of experiment [13] for (a) diatomic O2, (b) atomic 
Xe, (c) diatomic A"2, (d) atomic Ar irradiated by Ti:sapphire laser field of the wavelength 
A ~ 800 nm and Gaussian temporal profile (16) with the pulse duration r ^ 25 /s 
corresponding to about 9.43 optical cycles. The presented results correspond to various 
values of laser peak intensity /: (a) 2.0 • 10^^ W/cm'^ (dash-dotted line), 2.2 ■ 10^^ W/cm'^ 
(dashed line) and the saturation intensity 2.5 ■ 10^^ W/cm'^ (solid line), (b) 1.5 ■ 10^^ 
W/cm'^ (dash-dotted line), 1.6 ■ 10^^ W/cm'^ (dashed line) and the saturation intensity 
1.7 ■ 10^^ W/cm^ (sohd line), (c) 3.4 ■ 10^^ W/cm^ (dash-dotted line), 3.6 ■ 10^^ W/cm^ 
(dashed line) and the saturation intensity 3.75 ■ 10^^ W/cm'^ (solid line), (d) 4.0 ■ 10^^ 
W/cm'^ (dash-dotted line), 4.2 ■ 10^^ W/cm'^ (dashed line) and the saturation intensity 
4.4 ■ 10^^ W/cm^ (solid hue). 
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5. Fig. 5. (Color online). The differential harmonic emission rates Wj^ (k') [Eq.(13)] in 
laser-irradiated O2 and its atomic counterpart Xe vs the harmonic order A^ calculated 
under conditions of experiment [13] for the laser pulse of respective saturation peak inten- 
sity presently found and summarized in Table I. (a) The spectra of O2 corresponding to 
partial contributions from the Itt^ HOMO (open circles), the inner l7r„ (filled triangles) 
and Sag (filled squares) valence shells separately are presented vs overall O2 spectrum cor- 
responding to coherent contribution (32) from those three highest-lying molecular valence 
shells under consideration put together (stars), (b) The overall harmonic spectrum of O2 
presented vs the Xe spectrum (filled diamonds) calculated for respective (different from 
O2) saturation laser intensity. The results of Fig.Sb are to be compared with relevant 
ones presented in Fig. 5 of Ref.[13]. 

6. Fig. 6. (Color online). The differential harmonic emission rates wjy (k') [Eq.(13)] in 
strongly aligned N2 and its atomic counterpart Ar vs the harmonic order A^ calculated 
under conditions of experiment [13] for the laser pulse of respective saturation peak inten- 
sity presently found and summarized in Table I. (a) The spectra of N2 corresponding to 
partial contributions from the Sag HOMO (open circles), the inner Ivr^ (filled triangles) 
and 2cr„ (filled squares) valence shells separately and presented vs overall N2 spectrum 
corresponding to coherent contribution (31) from those three highest-lying molecular va- 
lence shells under consideration put together (stars), (b) The overall harmonic spectrum 
of N2 (stars) presented vs two harmonic spectra of Ar calculated for two different laser 
peak intensities - 4.4 ■ 10^^ W/cm'^ (filled diamonds) and 3.75 ■ 10^^ W/cm'^ (open squares) 
corresponding to the saturation intensities found for ionization oi Ar and N2, respectively 
The results of Fig.6b are to be compared with relevant ones presented in Fig.5 of Ref.[13]. 
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